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ON THE DIFFUSION OF HEAT IN A HOMOGENEOUS RECTAN- 
GULAR MASS, WITH SPECIAL REFERENCE TO BARS 
USED AS STANDARDS OF LENGTH. 

By Mr. R. S. WoopWARD, Washington, D. C. 

1. Considering the vast amount of time and labor which have been spent in 
developing the science and art of comparing standards of length, and the won- 
derful precision at present attainable in such work, it would seem somewhat 
surprising that more progress has not been made toward a complete theory of 
the behavior of metallic bars under varying temperature conditions. If we look 
over the scientific literature of the last half century, we shall find elaborate 
reports and memoirs on the appliances and methods adopted in determining 
coefficients of expansion and relative lengths of standards. But the chief object 
of all these appliances and methods has been to secure particular conditions, 
such as a steady low or high temperature of the bars under comparison; and the 
precise results for expansions and differences of length come out without refer- 
ence to the intermediate states and effects arising from the slow diffusion of heat 
through the bars while changing from one steady temperature to another. 
Indeed, a dexterous avoidance of all the more compiex phenomena presented 
by a cooling or heating bar has been the characteristic feature of experimental 
success. Our authorities on standards of length give us multitudes of facts and 
figures concerning the history of a bar during particular periods, but the history 
begins and ends abruptly and we are left to conjecture concerning the facts and 
figures of the more interesting periods. 

The reason. for this apparent anomaly, this paucity of information relative 
to the more recondite thermal properties of a bar, is obvious enough when we 
reflect that the present high grade of efficiency in the management of standards 
has only been attained after prolonged and diligent research. Obstacles quite 
independent of the standards themselves, of the most perplexing and baffling 
character, have had to be overcome. It is but natural and logical, therefore, that 
experimenters should have sought to avoid complexities as far as possible and 
direct their efforts so as to secure, in practical work particularly, the requisite 
precision at the minimum cost. 
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But precision is merely relative. What satisfies us to-day is almost sure to 
be unsatisfactory to-morrow. The perfection of methods and appliances which 
enabled us to overcome old difficulties has also revealed new ones, difficulties 
intimately related to, if not directly dependent on, the very states and conditions 
we have hitherto found it convenient to avoid or ignore. At present, precision 
of comparison appears to have reached a superior limit along the old lines of 
research; and future progress seems possible only by an advance in our knowl- 
edge of the more intricate thermal properties of metals. Apparently we must 
recur to the extraordinary studies of Fourier and Poisson; we must extend their 
work and adapt it to the practical needs of the computer; and we must learn 
how to determine readily for any standard the constants which define its intrinsic 
thermal relations and the relations it sustains to media in which it may be placed. 

2. It is proposed to consider in the following pages the problem of a rectan- 
gular bar cooling or heating in a medium whose temperature remains sensibly 
constant. This is probably the easiest of the problems presented in the work of 
comparing standards of length, but it appears to be also the one whose solution 
will be most useful in a practical way. 

Much, if not most, of the preliminary work has been done by Fourier * and 
Poisson*+, and the reader is referred to their treatises for all elementary details. 
In his Chapter VII Fourier has considered the case of a rectangular prism of 
infinite length, subjected to a constant temperature at one end, and exposed toa 
uniform current of air at zero temperature; and in his Chapter VIII he has given 
a pretty full treatment of the case of a cube cooling in air or any other medium 
of sensibly constant temperature. Poisson in his Chapter IX has discussed at 
great length the case of the distribution of heat in a bar whose transverse dimen- 
sions are very small, and in his Chapter XI he has briefly considered the case 
of a rectangular mass of unlimited dimensions. Neither of these writers, how- 
ever, appears to have had any practical applications in view, and elaborate and 
elegant as their work is, there remains much to be done to render it available for 
numerical applications. 

We shall assume that initially the bar has a uniform excess in temperature 
over that of the surrounding medium; that the conductivity and thermal capacity 
of the bar and the emissivity are constant within the range of temperature con- 
sidered. 

Let the length, breadth, and thickness of the bar be denoted by 2a, 26, and 
2c, respectively; and let any point within or on its bounding surface be defined 


* Théorie Analytigue de la Chaleur. Paris, 1822. Analytical Theory of Heat, by Joseph Fourier, 
translated, with notes, by Alexander Freeman, Cambridge, 1878. Our references in the text are to Free- 


man’s translation. 
+ Théorie Mathématique de la Chaleur. Paris, 1835. 
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in position by a system of rectangular co-ordinates, whose origin is at the center 
of the bar and whose axes of +, y, s are parallel respectively to a, 4,c. Call « the 
excess in temperature over that of the surrounding medium of any point +, y, 2 
at any time ¢ after the initial epoch. Then the partial differential equation, con- 
necting ¢, #, +, v, s, which obtains in all cases of the diffusion of heat, is * 


1 
(1) 


ct cs 


In this, & is the conductivity of the bar divided by its thermal capacity per unit 
of volume. 

In addition to the relation (1) pertaining to any element of the bar at any 
time ¢, there are three other equations which express the equalities between the 
heat conducted across any element of the faces of the bar and that carried away 
by the surrounding medium. If we call the conductivity A and the emissivity 
H, these relations are 


K 


u 
Hu for x - (2) 


My 9%) 


K~- = Ha for y= + (3) 


Finally, there is the initial relation expressing the fact that at the beginning 
of cooling all points of the bar have the same temperature. If 7, be the initial 
uniform excess in temperature of the bar above that of the surrounding medium, 


we have 
=U, for t=0. 


The analytical solution of the problem consists in finding a function 
u=F(u, x, ¥, 2, t), 


which will satisfy all of the relations (1) . . . (5). 

3. The complete integral of equation (1), for the present case, as shown by 
Fourier in his Chapter VIII, consists of the product of three infinite series and 
the initial excess of temperature #. We may express this integral thus: 


* See Fourier’s Chapter II, Theorem IV. 
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wherein 
A kth? cos As, 
v Bie COS 


"— cos vt. 


In these expressions A,, &,, C,, 4,, ft, ¥, are constants varying from term to 
term with the index ~, which is zero or any positive integer. These constants 
must be such that (6) will satisfy (2), . . . (5) as well as (1). 

lor brevity, put ,=-///K. Then, if we substitute w and its derivatives from 
(6) in (2), (3), and (4), there result 


ai, tan as, = ah =6, say, 
bm, tan by, — bh =, say, (7) 
cv, tan cv, =ch=£, Say. 


Since the values of ¢, 7, and £ are known, the equations just written will deter- 
mine the constants /,, 4, ¥,. Each of these equations, it will be observed, has 
an infinite number of real roots, the successive values of which are to be used in 
the several terms of 7,’, v,’’, v,’”’ in the order defined by the index x. 


The constants 4, 4,, C, are determined in the following manner: Equation 
(5) requires that when ¢=0, 


that is, considering v,’ alone, 
1 = A, cos Ayr + A, cos Ax + A, cos Ax + 


If we multiply both sides of this equation by cosd,rdx and integrate 
between the limits o and a, it will be found that * 
sin ad, 

A, = 2a4,, + sin 2aA, (8) 

Corresponding expressions for 2, and C, may be written down by an obvious 

interchange of symbols. 

4. It is seen, therefore, that the use of equations (6) as they stand, requires 

the roots (or a certain number of them) of the transcendental equations (7), these 

roots entering respectively into each of the three factors of the corresponding 


*See Fourier, Chapter VII, p. 313. 
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term of v,’, v,/’, or v,’".. Herein lies the chief obstacle to the application of (6), 
and the principal feature of our work consists in showing how this obstacle may 
be overcome. 

We shall show first how the roots of (7) may be found, though our solution 
of the problem requires in any case only the smallest of these roots. 

For brevity, let us write 


Then (7) become 


It will suffice to consider the first of these equations, since they are all of the 
same type. 

§ = ah is of zero dimension as respects length*, and it may, in general, have 
any value between 0 and o. 

In the case of an iron bar cooling in air, 2 may be about + with the cen- 


timeter as unit of length; so that for a bar of iron one meter long, we should 


have approximately (roughly, perhaps) ¢—,'y. But for other media, 4 may be 
so large as to make ¢ exceed unity when a is no greater than a half meter. In 
any event, ¢ will exceed unity if a is sufficiently great. It becomes essential, 
therefore, to distinguish two cases; namely, that in which ¢<1 and that in 


which ¢>1. 


First, consider the case wherein $<1, with a view to expressing @, in a 
series of ascending powers of ¢. When ¢ =o}, we must have 0), — uz; i.e. 0, is 
either zero or a multiple of z. Differentiating the first of (10), there results 


dat, 


as 


This gives, for § =o, 
1en #7 =O, 
(nz)! when 
* See Fourier, Chapter II, p. 128. 
+Sir Wm. Thomson gives (Encyclopedia Britannica, gth Edition, Article Heat) 4 — 0.0002 for a 
polished copper surface, and A= 0.16 for iron. Assuming the above value of 7/ to apply for iron as well 


as copper we have 4 = ;}5. 
} The physical meaning of this condition is that there can be no escape of heat from the bar. 
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ee 
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tan ¢, = 4, (10) 
tan ¢, = &. 
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It appears, therefore, that we cannot apply Maclaurin’s series to the development 


of the first value of #,. But if we substitute for tan 4, its value in series, we have 


whence, by reversion, we get for the first roots of (10) 


2 1 22 4 <3 16 =4 JT * 
> 3” 45” 9 45 » Ts. 
| 4 v3 _ 16 «At 
Co 34 45% 945% (12) 
f, 2 1 #2 4 #3 


These series will suffice for the calculation of the first values of 4,, ¢,, ¢,. The 
convergence is rapid for small values of ¢, 7, £, and quite adequate even when 
they are as great as unity. 

For the larger values of 4, we may apply Maclaurin’s series. Thus, by 
eliminating the trigonometric forms from (11), we get 


as 


t 

— 


whence, by an application of the theorem of Leibnitz, we readily find that, 


for | =0, 
d= if ) 
a=? 2( 777) (13) 
9( 
d= 2(uz)”, 
Hence, we write 


n—1, 


This series converges rapidly, especially for large values of 2; it shows, in ‘act, 
that the larger the root the more nearly is it equal to zz, 
Second, when ¢>1, we may write the first of (10) in the form 
cot 6. 
wherein (15) 
Now the first of these equations gives for ¢,—=0,t 


* We shall use throughout this paper 7;, to indicate that the next term of the series beyond those 
given is aterm of the mth order. Thus the next term in the first of (12) is of the 5th order with respect 


to the variable &. 
+ The physical meaning of this is that the heat is instantly dissipated on reaching the surface of the 


bar, or escapes “ freely.”’ 
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1)z7= J, say, 


Therefore, in accordance with Maclaurin’s series, 
6, = [1 + — + J). (17) 
+ 


Except for large values of =, this series converges too slowly to be of much 
practical value. Our process, however, does not require its use. We give it only 
to facilitate a numerical verification cf the equivalence of the original and trans- 
formed expressions for 1’, ’ 

5. It will appear in the sequ-l that, for small values of §, 7, and €, the series 
of (6) converge rapidly toward their first terms with the lapse of time from the 


See § 7. 


initial epoch. On this account, it is desirable to have a ready method of com- 
puting 7%’, %’’, and v,/’’.. The factors in these terms are all in convenient form 
except d,, 4), and C,, which are implicit functions of , 7, and £, respectively, as 


defined by (8), (g), and (10). They may be expressed directly in terms of ¢, 7, ¢ 
in the following manner. Taking A, as the type, (8) and (9) give 
4 sin 4, 


sin 20, 


A 


Substituting for sin #, and sin 2, their values in series, we find 


‘ 4 
a)* cos 4, 
a 


= By — cos 


as; 
as, 
as, 
: 
| 
* 
| 192 1 : 
Replacing in this by its value i here resul 
eplacing (/,° in this by its value in (12), there result a ae 
7, 
From (6) and (g), we have 
(19) 
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Computation of numerical values from these expressions is now easy, since 
Yo. are given by (12) and A,, CG by (18). 
In addition to 7,’, 7,’’, and 7,/’’ we need 


wy 


~ Uns On» 


1 nN H 


These we proceed to expand as functions of ¢, 4, and £, respectively, by Mac- 


laurin’s series. 


By reference to (8), (g), and (14) it is seen that when »>1 and ¢=0, 
A,, O, 
e Ath? At(nw a)? 


5 
COS = COS ; 


a 
aA, 
as 
fA. 
3° = 
as 
di, 
as 
Hence, when = = 0, 
( 3 ) e COS NI, 
as a 
n+l 
12(—1 x 
a= a 
8h ir x 
ann a 


n+l 
41 I A 

( ) e kt(nw a)? sin mz, 
QT HN a 


If now we write for brevity, 


nm n+l 
= ( ) e kt(nn /a)? 
nH a 
2 I n+l x (20) 
Tv nN a 
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we have tor =o, 


Finally, if we denote by S,”’, S,’’, etc., and by S,’’’, S,’’’, etc., the series 
y, 1 2 1 2 

corresponding to (20) when a and + are replaced by @ and y, and ¢ and sz, 
respectively, there result 


35,/ 


T 


Equations (12), and (18) to (21) afford a complete solution of the problem for 
the case wherein ¢, 7, £ are each less than unity. The series in (20) converge 
with such rapidity that for any but the earliest stages of cooling their first terms 
will suffice. For practical applications to a standard of length of ordinary 
dimensions, cooling or heating in air, we may neglect all terms after the first 
in (21). 

The values of 5S,’, etc., in (20) might be expressed by a series of definite 
integrals similar to that derived in §8, but it does not appear that any practical 
advantage would be gained thereby. 

6. Considering still the case in which , 7, and £ are each less than unity, it 
will be of interest to have expressions for certain average temperatures, namely : 
(a) the average temperature of the whole mass of the bar, (,4) the average tem- 
perature of any face of the bar, and (j7) the average temperature along its longi- 
tudinal axis. 

(a). The average temperature of the whole mass is defined by the definite 
integral 

+a +e 
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A glance at (20) and (21) shows that this is equal to 
a +6 re 


a 


whence by reference to (19) we get the following result: 


_ sin sin ¢, sin 


(a)=4,A,B,C 2 


! 
which, by reason of (12) and (18), is easily reduced to 


(4). The average temperature of any face of the bar, for example the face 


2a « 26, is expressed thus: 
a 
. 


a 
wherein ¢,/’ stands for the second member of the last of equations (21) when 
c. Integrating, there results 
( ale Me?) 7,). (23) 


It is easily seen that this differs from (22) by a term of the first order in &. 
(7). In a similar manner we find for the average temperature along the lon- 
gitudinal axis of the bar, 


(7) —= My 73), (24) 


in which «,’’ and a,’’’ stand respectively for the second and third of (21) when 


y==oand s=o. Obviously, (24) differs from (22) and (23) by terms of the first 


order ir 7% and £. 
7. Returning now to the case in which ¢, %, and £ are each greater than 
unity, we proceed to show how the sums in (6) may be expanded in series of 


ascending powers of $~',7~', and £~'. Referring to (15) and (16) and using 


for brevity 4, ,3, 7 as defined below, we get for |~'=¢,=0, 


a= A, = + 2(—1)'9-}, 
ou 
ca n 
>] 


* Compare Fourier, Chapter VIII, p. 328, for the case in which a=6=c. 
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a)? 


2ht (9 


[ a) - (240)*(a)* Je 


y=cos cos 
a 

x 

sin 

a a 


— 20— sin 
a a 


If now we put I,’ for when the preceding relations enable 
us to derive readily the following: 

en 


At(2m (w ~2a)* cos l(2n j 
I a 


Hence by Maclaurin’s series 


pu — 


in which l,’’, V,/”" and the derivatives of V”’, V’”’ with respect to 7~' and £7! 
are obtained from (25) by replacing a, x by 4, y and c, z respectively. 

The applicability of (26) for numerical computations will depend chiefly on 
the rapidity of convergence of the series expressing Il’,’, V,/’, and V,/’".. These 
series (see the first of (25) ) converge rapidly when the time is sufficiently great, 
but rather slowly during the earlier stages of cooling. They may be expressed, 
however, by a series of definite integrals which converge with extreme sud- 
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denness for all but large values of the time; and as the process of transformation 
from the former to the latter series is of great utility in other problems as well as 
in this, we shall give the steps of its present application rather fully in the next 


section. 
8. If in the equation * 
3 20. 3.3 
| 
we put 
2 
a = 
and 
the value of V’,’ in (25) becomes 
3 cos 4(2” +1)z-—-cos 4(2"-+ 1)z 
0 ah / 2 
(747), 2n a 
1 ( 
cos 4(2a-+1)z | 
2 ( 1) 2 a (27) 
my 


Now the series under the sign + in this equation is periodically equal to zero 
and 37, as 7 varies over its range from 0 to «; and hence IV,’ consists of a defi- 
nite integral whose successive parts are multiplied alternately by o and 4z. 
These facts may be shown in the following manner: 
Take the difference between the well known equation 
(1 


cos CSP +s cos cos 3p 4 


and the equation which follows from the substitution of z+ f for ». The result is 


i—z* 
} I1—2¢ cos I 22 cos p § Cos p+ COS 3P + 


COS P + 


For in this, substitute go° —(G We find 


2g sin(G + p)+ 1+ 2¢ sin (G + 
g sin (G+ p)—g' sin (G +p) + 


* See Vol. III, p. 78, of the ANNALS OF MATHEMATICS, 
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Multiply this last by dp and integrate between the limits 
result is 


n+P 


I 2¢ sin (G Si V, say. 


P 


But Q is evidently identical with the quantity under the sign ¥ in (27) when 


2a 


The element-functions of the integrals in the first member of (28) are obvi- 
ously zero for g=1 except when G + f is $(4/ + 1)z in the first and }(47 + 3) = 
in the second, 7 being zero or any integer. In these exceptional cases the ele- 
ment-functions are infinite; and the values of the integrals depend on the limits 

~Pand P += only so far as they determine within what range 7 can give to 
G+/p the above named multiples of }z. Hence, in evaluating either of the 
integrals, as for example the first, for any case in which sin (G+ f) = + 1, we 
may enlarge the limits to —« and - © without affecting the result. This 
premised, the required evaluation may be accomplished by a process frequently 
employed by Poisson in the treatment of similar integrals.* 

Substitute for g, 1 —w, where w is a small quantity approximating zero as 
g approaches unity. Also, replace G + p by $(42 + 1)z + @ in the first integral 
and by }(42 + 3)z +g in the second, g being, like w, a small quantity. Then, 
neglecting terms of the third and higher orders in g, either of the integrals 
becomes 


d 

Wag 
q 


2 


wherein g, and g,, as seen above, may have any values greater than an infini- 
tesimal. Hence the value of either integral is 


72 
7 —arc tan 
w 


4 [ arc tan 


when 


* See, for instance, Chapters VII and VIII of the 7héorte Mathématique de la Chaleur. 
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We must now find the limits of G, or rather 7 as defined by the second of 
(2g), within which the factor Q is }z. For the first integral the lower limits of 
7 are given by the equation 


C= 1 (42 + 1)% 


or y= (42+ 1)a— 
whence 

y== —( a+-2x) forz=o, 

- (3a—-) forz=1, 

=—-+(7a—-+) for 7=2, 


and the upper limits are given by the equation 


whence 
for z=o0, 
for 7=1, 
=oa-+ for 7=2, 


Similarly, the lower and upper limits for the second integral are, respectively, 
a—-x and 3a-- « for 7=—0, 
and 7a for7=1, 
Hence, bearing in mind that 7 is positive only, we find for the transformed 


expression of (27) 


aa 3a +2 sa +2 
o a-—x 3a — x 
(30) 
a—x 3a—x 5a— x 
a x 
If now we write for brevity 
and replace 7 by 27, (472), there results 
3%0 — 7 
it = -| dy .. (32) 
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The corresponding expressions for ’,'’ and V’,'”" are obtained by merely inter- 
changing a, with 4, y and ¢, 

It will be observed, of course, that the series of definite integrals in (30) or 
(32) possesses all the properties and requirements of the less complex series in 
the first of (25). We shall not trace out these details, but the reader’s attention 
may be drawn to the remarkable way in which (30) and (32) become zero for 
4= and unity for ‘=o. 

In case ¢, 7, € are so large t':at we may neglect all terms after the first 
in (26), 


This is always an approximate expression for the temperature when §¢, 7, 


and £ are greater than unity; and, other relations remaining constant, its degree 


of approximation increases as the dimensions a, 4, ¢c increase. For points near 


the surface of a very large mass, comparable in dimerisions and other circum- 
stances with the earth, say, we may write 
(a (6 —y) ay (at) ( ay (kt 

u = 8u,(z)—! e Ydy e 

° 

For points near the surface of a slab of infinite dimensions in the directions 
of 6 and c, and finite but great dimension in the direction a, the last equation 
gives 
(a— x) (At) 
fe 


7% 


a well known result. 
The derivatives of l’,’ with respect to x and ¢ which enter the last two 
equations (25), are easily found from (32). For brevity, write 


- 2ur, 


= 
f 
4 R, a f I ) [< Py 
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Then we have 


on 
R, 
Rey 
R, 
ct ty az 
I 
I 
at) — 244). 


ut 
— 


g. We are prepared now to derive any species of average value of the tem- 
perature for the case under consideration. Of these, two may suffice, namely: 
(a) the average temperature of the whole mass and (,4) the average temperature 


over any face of the mass as 2a x 20. 
(4). For brevity put 


=(2a) 


a 


M” =(26)-'{ "Vay, 


(2c)-* 


Then 


If now we write 


it appears from (25) and (26) that ~ 
7, 


Similar expressions obtain, of course, for 4/’’ and 1/’”’. 


(36) 


(38) 
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It remains to evaluate the integrals in (37). They are simple in form but 
present some of the difficulties peculiar to the class of functions to which they 
belong. The first of them is readily evaluated by direct operation on the first 
of (25), while the second comes by a more tedious process through the inter- 
vention of the expressions in (31), (34), and (35). The results are 


kt(2n + 4a? 


(39) 


(2% +1) 


The steps in the derivation of , are briefly these: 


vo 


2) (At COR tre 


ay 
vo 


+-1)% 


(2ur, + — p)dp 
e 


1)%o 


whence the expression in the text. 

The second members of (39) though complex in form are very easily evalu- 
ated, by reason of the rapid convergence of the series on which they depend. 
It may be noted that for ¢= 0, V—1 and .\,=—0, and that for ¢-—-», NV 
N, =0, as required by the conditions of the problem. 

(3). To get the average temperature of the face 2a « 26, we have only to 
multiply ~,J7’/" by the value of I’’” in (26) when z==c. Calling this value 
we have 


c 


(3) V2" (40) 


10. It is believed that the solutions given by (6), (21), and (26) will answer 
most of the requirements in practical applications to standards of length of rect- 
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angular form. It is possible, however, that the developments with respect to 
§,7, and £ may prove to be in some cases insufficiently convergent. This can 
happen only when ¢, 7, and £ are near unity, and the obvious modification requi- 
site for the treatment of such cases, is a development with respect to 1 —&, 1 —4, 
and 1—%. But in the absence of a definite knowledge of the circumstances 
giving rise to this phase of the problem, it does not appear useful to do more 
than suggest the modification. 

11. For the purpose of illustrating the application of the preceding theory, 
we shall work out in detail some numerical examples. 

First, take the case of an iron bar, 1” 100” long and 4” X 4” in cross 
section, cooling in air. We have for iron, according to Sir W. Thomson,* in 
C. G.S. units, #- -0.185, and we may, as in §4, assume 4=-%},5. Then, since 


90", 4-2", 


1 


1 
is % T00° 


With these data, equations (12) give 


log 4° - - 8.78689 — 10, 


log ¢,° 7.39758 — 10, 


log 7.39758 — 
Also 
10’ 
2 5! 49, 
2 51 49; 
and 


log & 8 4.65621 — 10, 

| fo | =: 606269 — 10, 

log & 6.06269 — 10. 


From (18) we find 
log A, = 0.00440, 
log B, = 0.00018, 
log C, = 0.00018. 


Encyclopedia Britannica, 9th Ed., article Heat. 
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We are now prepared to compute the factors 1”, 1’, and 1” of (21) for 
any values of ¢, x,y,z. For example, let us determine these factors and hence 
the temperatures for several points along the axis of the bar at the end of twenty 
minutes from the initial epoch. That is, make 1200"; successively 0, 
ha, 3a, and a; and y= z= 0. From (19) and (20) we get by means of the 
preceding data the following table of values for use in (21): 


1.0047 + 0.0828 0255 5 + 0.0000 


1.0028 | + 0.0596 +oo18 53.) +00170 | 
0.9970 | 4 0.0015 0.0001 | + 0,0268 
0.9875 | —0.0596 | - 53 | 0.0208 
0.9741 | —0.0859 —0.0257 | 75 | 0.0000 


In addition we have 
0.87054, 


0.0828, 


and 


Hence we find the following table of values of I’, in which the terms in ¢ and = 
are given separately in order to show their magnitudes relatively to 7,’ : 


= 0.9996 

+ 0.00016 0.9993 
0.00010 0.9970 
0.00001 0.9912 
| 


Blo 


| 
‘ 
or? 2 0.03 30. : 
a | > As |’ 
Oo, + 1.00 
+ 1.0028 | — 0.0037 | + 
| + 0.9970 | eae 
+ 0.9741 | + 0.0054 | — 
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We find also that 


— = 0.87054 — 0.000207 + 0.0000002 


0.87033. 


It is now only necessary to form the products ’1”’V’”’ in order to get the 
temperature w# in terms of #,. The results are 


a u 

a My 

0.7571 

1 

4 0.7 579 

0.7553 

0.7509 

I 0.7420 


It appears from the values in the second paragraph above that the terms in 
=, 77, and £? are, for the assumed values of those ratios, already practically insig- 
nificant at the end of twenty minutes from the initial epoch. Hence in computing 
temperatures for greater values of the time we may properly neglect all terms of 
the second order. For the purpose of showing the march of the temperature 
with the time, we have computed the values which w has at the above named 
points along the axis of the bar for intervals of twenty minutes up to ¢= 5” 40”. 
Since minute precision is not needed in this work, we have used ” = v,' — 5,’ 
and |”’= |”"’ —7,/’.. The resulting values are shown in the table below. The 
average temperature of the whole mass of the bar obtained from (22) is also 
given in the last column of the table. This average is only slightly greater 


(;slypth part) than the average temperature of the face 2@ « 26.* 


* The question whether under similar circumstances there can be any detrimental internal strain due 
to inequalities in distribution of temperature in standards of length of ordinary forms and dimensions, 
appears to receive a decisively negative answer in the computed values of the table, Such standards are 
generally smaller in cross section than our assumed bar, and would, therefore, other things being equal, 
sooner attain a sensibly uniform temperature. Hence we infer thet, so far as flexure due to temperature is 
concerned, it is immaterial where the graduations on a standard are placed—whether on its external surface 
or on its neutral plane, 
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Values of temperature # on axis of bar where x, a Average | 
Time for | 
t mass of 
ra) I bar. 


A m 
O 20 0.757% O.757M, | 0.751%, 0.742%, 0.754%, 
O 40 0.573 0.572 0.570 0.565 0.558 0.568 
I OO 0.433 0.432 0.430 0.426 0.420 0.428 
I 20 0.326 0.326 0.324 0.321 0.317 0.323 
| 1 40 0.246 0.246 0.244 0.242 0.239 0.243 
| 2 00 0.185 0.185 0.184 | 0.182 0.180 0.183 
| 2 20. 0.140 0.140 0.139 0.137 0.136 0.138 
2 40 0.105 0.105 0.104 0.103 0.102 0.104 
| 3 00 0.079 0.079 0.079 0.078 0.077 0.079 
| 3 20 0060 0060 0059 | 0.059 0.058 0.059 
| 3 40 0.045 0.045 0.045 0.044 0.044 0.045 
| 4 00 | 0.034 0.034 0.034 0.033 0.033 0.034 
| 4 20 0.026 0.026 0.025 | 0.025 0.025 0.025 
4 40. OOIg 0.019 0.019 | 0.019 0.019 0.019 
5 00 0.014 0.014 0.014 0.014 0.014 
5 20 0.011 | 0.011 0.011 
5 40 0.008 0.008 0.008 | 0,008 0.008 0.008 


12. Turning now to the question of numerical applications for the case in 
which ¢, 7, and £ exceed unity, it will be advantageous to first set down some 
additional formulas for computing the differential coefficients which appear in 
the second members of (25) and (26). 

To avoid writing frequently recurring complex expressions, let us make the 
following abbreviations: 


nN 
(—1)" (2% + 470)" 2" Cos + 1) za~' 4, 
nN Oo 
X, (—1)* 7470)" +1)” sin + 1) 2, 
N= (41) 
As (—1)" (2a + +1)? cos + 1) 4, 
nN 
n 
(—1)" (24 + sin (2% + 1) 2, 


ro? =a? ght. 


; 
‘= 
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Then, the second and third of (25) become 


Lar, * X, + 2a-'"2X,, 

(42) 
2.2 1 3, -2) at, 


The equations just given will be useful when 7, is small, since the functions 
in (41) will then converge very rapidly. For the converse case, we must make 
use of the differential coefficients given by (35). Introducing these in the second 


and third of (25), we get 


(rk, + R,) 
2) ad vr) 
4 + 7°) R, + ark, + R,). 
esp 


r? 


4ht 
The last expression assumes very easily managed forms for certain values of r 
used in the sequel. For the benefit of the reader who may desire to verify our 


work, these forms are given below: 


a 
C3) 
4 3 25¢ -) for r= 4%, (44) 
0 47% 5? 1675? | 6% 5? 
04 (« ge 3 ret, 


The series in (43) and (44) converge with extreme suddenness for large 
values of 7, and they cease to be applicable only when 7, is very small. 

With preliminaries thus arranged, let us consider the behavior of a bar of 
the same dimensions and material as the bar of the previous example, cooling in 
a medium which produces an emissivity great enough to make §, 7, and £ exceed 
unity. Without assigning to these symbols any particular values, it appears that 
the bar would very quickly assume the temperature of the medium. For we 
have, as in §11, 

50", = 2", 


and k= 0.185; 
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and since for the transverse section of the bar 


kt 


the quantities and the derivatives of 1”’ and relatively to and 
£~' become practically wz/ soon after the initial epoch. Thus, from (25) and (42), 
we find the following table of values which show the rapid decrease of those 


quantities with the time 


0.041 + 0.991 
0.001 3 0.018 0.195 


| QO 0.00004 + 0.0009 | 0.016 


Hence we see that, with such values of =, 7, and £ as render equations (26) appli- 
cable, the bar in question would have sensibly the same temperature as the 
surrounding medium in less than two minutes after the initial epoch. 

Other conditions remaining constant, the rate of cooling of a mass dimin- 
ishes very rapidly with an increase of its dimensions. In order to strikingly 
illustrate this fact by means of the analysis, we may contrast the cooling of the 
bar just considered with the cooling of a cubic meter of the same material. The 
use of the formulas will be sufficiently indicated by computing the temperatures, 
corresponding to a few values of the time, at the five points designated by the 
following values of their co-ordinates: 


To compute V’, V’”’, ’”’ for such values of the time as we have assumed, 
one may use equations (25) and (42) or equations (32) and (43), the latter being 
best for the smaller and the former best for the greater times. We have used 
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both sets of formulas as a check on the computed results. The quantities needed 
in equations (26) are given in the following table. V,’, V,/’, and V,/” and the 
corresponding derivatives of V’, V’’, and V’” are equal for += y= s=0. 
Hence the values in the first section of the table are those required by the 


second and third of (26). 


| 
A= AE 3 V,! AEs 
O 20 | 0.965 0.227 2.55 0.764 0.940 5-33 
| 00 | 0.659 0.852 3.50 0.467 0.874 - 1.85 
2 0.342 0.905 4.13 0.242 0.825 | 1.11 
3 0.177 | 0.699 1.26 0.125 0.587 | 3.80 
a, y=0, s=0 
| 
ye OV 
0 
Cy] 
| h | 
| O 20 oO 1.893 0.00 
| 1 OO O 1.042 0.84 
| 2 00 O 0.537 + 5.45 
| 3 00 O 0.278 9.37 


Let us now suppose that 74 10; te, 10, or h 
H/K--+. Then we have for the centre of the cube when ¢. = 20”, 
— = 0.965 4+- 0.023 — 0.013 + 7; 
0.975, 
and U (0.975 == 0.927 Uy. 


For a point on either axis of the cube half way between the centre and a face, we 


have when 1", 


Vy’ 0.467 + 0.087 — 0.009 + 73, 
— V'" = 0.659 + 0.085 — 0.018 73; 
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and hence 


== 0.545 X (0.726)?u, = 0.287 


In a similar manner the temperature at the centre of a face of the cube is 
found. The following table gives all of the temperatures derivable from the data 
in the preceding table. The values of the co-ordinates may, of course, be inter- 
changed, so that the temperatures given in the last two columns are common 
respectively to six points. 


Time. Temperature. 


At middle 
At centre of At centre of 
”) 


4 face 
sem 


cube. 


0.9274, 0.180%, 
0.384 0.288 0.057 
0.069 0.054 0.014 


0.016 0.013 0.005 


As a final example in this section we may give the average values of the 
temperature of the whole mass and of any face of the cube just considered. 
These values are easily computed from cquations (36) to (40). Using the same 
values of the time as those in the last table above, we arrive at the following 
results, exhibited also in tabular form. The quantities V and .V, are computed 
from (39), 17’ from (38), and V’”’ comes from the data in the third section of the 


second table above. It is assumed that § = 7 = = 10. 


Of any face. 


0.689 | + 1.34 0.56%, 0.13%, 
00 | 0.422 | + 1.53 UT; 0.19 | 0.04 
00 | 0.331 | + 1.37 | 0. 0.10 0.02 
0.113 | + 0.99 0.01 0.01 
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13. It appears desirable to indicate in the last section of this paper in what 
ways the conditions of the problem considered and the analysis derived may 
be useful in determining the thermal constants of a bar; namely, the thermal 
capacity of C, the conductivity A, and the emissivity H. 

The principal conditions of the problem we have treated are first, that the 
bar has initially a uniform temperature, and second, that it cools or heats in a 
medium whose temperature remains sensibly constant. The facts most readily 
observed under these conditions ‘are the initial excess in temperature of the bar 
above that of the medium and the variation of its temperature or some function 
thereof with the time.- From a sufficient number of such facts it would appear 
that all of the constants might be determined. To fix the ideas, suppose the 
excess in temperature at the centre of the bar or the average excess of its entire 
mass, can be measured at any time. Call this excess vw’. Then we have 


u'’=/(C, H, K); 


and if we write 


=C, 
H,+ dH=H, 
K,+ I4K=K, 


in which C,, ,, and A, are approximate values of C, H, and X found in any 
manner, we shall have the following observation-equation: 


IC TOK tT (Cy, Hy, Ko) — = Aw’. (45) 

From a series of such equations, differing from each other sufficiently, one 
could by the .method of least squares derive precise corrections to the approxi- 
mate values C,, /7,,and A,. The applicability of this process, however, is limited 
narrowly by the difficulty of securing a sufficient range in temperature in a given 
bar, or by the difficulty of securing a sufficient range of dimensions in different 
bars or masses of a given material. 

Although it does not appear practicable under the special conditions of the 
present problem to derive all the thermal constants of a bar, there exist between 
those constants under certain circumstances some relations which may be accu- 
rately determined. ‘These we proceed to point out. 

Firstly, when the emissivity and the dimensions of the bar are smali, we 
have from (g) and (22), the average temperature of the bar’s mass 


uu! = ue — + 0%) + 


whence + + log e = log 
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and from (g) and (12) 


42+ + +67! 


Hence, to terms of the first order in we have 


(a~' + 6-'+7') log e= log (46) 
This expression gives the emissivity 47 in terms of the thermal capacity C when 
the second member is known by observation. Under the assumed condition 
that 47 is small, w’ may be safely inferred from the length of the bar (supposing 
it to be of small cross section). In case tnere are several observed values of x’ 
and the corresponding time ¢, let W be the most probable value of the second 
member of (46). Then, since #’ or log w’ only can properiy be considered sub- 
ject to error of observation, we have in the usual notation of least squares 


u 
E log 


ay" 


From this and (46) inerefore, 


CW(a- + 4 


log e 


(48) 
Secondly, when the emissivity is large or when ¢, 7 and £ are such, we 

have from (6) and (26), 

Now V,/, V,'’, V,/" and the derivatives of V’, V’’, V’” involve the coefficient of 
diffusion or C and K only; and hence when these are knowa an observation 
of « and ¢ will give // through the relation (49). On the other hand, if ¢, 7, and 
¢ are so large that the terms beyond the first in (4g) may be neglected, an obser- 
vation of « and ¢ will give the coefficient of diffusion &. The method of treating 
a series of observed values of uw and ¢ to determine # or 4, is sufficiently indicated 
by the remarks pertaining to equations (45) and (47). 
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NOTE ON THE THEORY OF IMAGES. 
By Dr. Ro_itin A. HARRIs, Jamestown, N. Y. 


Article 36, found upon pages 84-6, Vol. TV of this journal, is susceptible of 
the following interpretation :— 

In that article nothing was said concerning the properties of and 7. It was 
implicitly assumed that, like ordinary algebraic quantities, they were commu- 
tative, associative, and distributive; and that 7? = 7? = — 1. 
In expanding ¢ wy ++ 7z) we have 


X +4V¥ + 


Now if we confine ourselves to space of three dimensions instead of four, we 
project the latter into the former by suppréssing the term 7W; i. e. by omitting 
all terms of the form Ki"s"y"s" where m and n are both odd. , 

By this interpretation what was called the path of U is really a projection of 
the same. 

If 7, 7, (= 7) be commutative, associative, and distributive, then generally 
g(x + wy -+ sz + kw) is developable in powers of zy + 7z + Aw and is 


7Z+- kW; 
and in four-dimentional space the locus of the point xyzw has an image, the 


locus of XYZ. Now if we suppress &W we project the path of U, the image 
of the path of + + yz + 0), into the (1, 2, 7)-space. 


Professor Oliver suggests that by writing a=1, 1,7 S=JV —1, 
0 = —&, this (1, 2, 7, 4)-algebra takes the symmetric form 


as fa = OF, 
ay = fa = = po. 
ad = da = fy = 78. 


3 
; 
J 
| 
] 
} 
4 
2 
4 
i 
; 
ur 
| 
q 
4 


KUMMELL. FUNDAMENTAL THEOREMS OF MENSURATION. 129 


ON SOME FUNDAMENTAL THEOREMS OF MENSURATION IN 
ONE, TWO, AND THREE DIMENSIONS.* 


By Mr. C. H. KUMMELL, Washington, D. C. 


This article has for its subject the following three analogies: 


~ 


The first of these formule gives the distance of two points whose distances 
from the origin +, and x, are given. 

The second gives the area of any quadrilateral in terms of the co-ordinates 
of its vertices. 

Similarly, the third gives the volume of any octahedron, or any solid within 
six points in space, in terms of orthogonal co-ordinates. 

It will be noticed that in place of the usual fundamental forms, distance, 
triangle, and tetrahedron, we have here distance, quadrilateral, and octahedron, 
and it will be easily seen that above general formule can be used for triangle 
and tetrahedron by causing one or two points respectively to coincide. This 
does not, however, essentially simplify the formule. We have thus for the area 
of triangle 1. 2. 3 


42 — 4, 4%, — 


(4) 


Is 


and a variety of other forms obtained by transformation. 
We have also for the volume of the tetrahedron 1. 2. 3. 4 


== £60. 


* Read before the Mathematical Section of the Philosophical Society of Washington, D. C. 
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On the other hand, it will appear from the demonstration of these formule 
that neither the area of any other polygon nor the volume of any other polyhe- 
dron can be expressed by a single determinant. The area of any polygon of x 
sides is the sum of # trapezoids formed by each side with its projection on one 
of the axes, taking it positive if in following the contour we go farther from the 
other axis, and negative if we approach it. The area of the trapezoid on the 
side 1. 2 is = } (4, — %) (32 +), and that of the other trapezoids is obtained 
with the proper sign by simply circulating subscripts in the order: 1, 2, 3, ..- . 
n,1. Thus we have 


A= (42 — (2 +1) + — 42) + 2) + — 45) (Me + Is) 
Vor V3, Vy Is Vs, In 
2 2 


If ~ is even, the area is thus expressed by 4 determinants, which are not com- 
binable, in general, because they have no common column. In the case of the 
quadrilateral, § coincides with 1, and we have 

4) — 43, 


Fe Fe Ji — V3 


In using this formula care should be taken that the elements of this determinant 
are the projections of ¢zterna/ diagonals or in case of a re-entrant qualrilateral at 
least of one internal diagonal, for if we take the points in a different order we 
obtain the area of a crossed quadrilateral. If these rules are attended to we 
obtain the area of a quadrilateral in the ordinary sense except with regard to 
sign. If we wish. it positive we must figure the vertices so that the positive 
diagonal product (#, — +,) (7, — J) is positive and as large as possible. 

Formula (2) is not essentially new, except, perhaps, being expressed by a 
determinant, for it may be easily recognized ina rule for computing area as used 
in land surveying.* 

Formula (3) may be demonstrated in an analogous manner. For any poly- 
hedron may be expressed by the algebraic sum of as many trapezoidal prisms 


* See, for instance, Davies and Peck’s “Dictionary of Mathematics,” Article Survey, p. 550. 
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as it has faces, taking those negative which are wholly external. It is of course 
immaterial on which co-ordinate plane the faces are projected, but we shall use 
the horizontal plane. 

Suppose we wish to have the volume of a trapezoidal prism formed by pro- 
jecting the face (m, , 0) on the horizontal plane. It is well known to be = 
projection < mean height, or, using formula (4), we have 


or we may write 


3! (6) 


where the third row may also be written in any order whatever, because their 
minors are equivalent. If we reverse the order of the points, we have 


EA ny Amy 4m — 


I 


Now it is clear that if a trapezoidal prism is external to a polyhedron, then 
the order of the points on the latter is the reverse to that on the former. If, 
therefore, in using (6) for the trapezoidal prisms on each face of a polyhedron 
we take the points always in the order as they appear on the latter, we obtain 
their volumes with the proper sign, and their sum total will be the volume of the 
polyhedron. 

To aid the perception I have constructed a stereoscopic diagram of an octa- 
hedron, exhibiting also the eight trapezoidal prisms Piss, 
Pris. Ps 6.5 Pa 2.6, and we have the volume, 
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I 
~ ~ 
“1 “6 
Be 45, — 43, — 4%, 
I | 
V4 Vs — Vn Va — Ia 
“39 
I 
V4 3! Ja — Ia V2—- Ie Ve— Ie 
Ze 


Arranging this in terms of 2), 22, 25, 24, %;, and 2, we have 
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If the co-ordinates are measured along oblique axes, so that the y-axis inclines 
to the z-axis by the angle y,, and the z-axis to the +y-plane by the angle z,,, 


then, if we replace in all preceding formule 
y by y sin y,, 


by 2 sin 2,,; 
we have the formule for oblique co-ordinates 
I | 
1 |%— . 
A=— sin y,, 
25 I3— Si | 
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and if the a-axis coincides with diagonals 2. 4 of the quadrilateral, and 3.6 of 
the octahedron, the y-axis with diagonals 1. 3 and 5. 2, and the s-axis with diag- 
onal 1. 4, then the elements in the positive diagonal of the determinants are the 
respective diagonals, while all side elements vanish, and we have 


L= (1. 2), 
1! 

A= 4) (1. 3) sin y,, is") 

V= ~ (3.6) (2. 5) 4) sin y, sin 2,,. (3) 
a* 
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Snow that the locus of the centres of equilateral hyperbolas circumscribed 
to a given triangle is the nine-points circle of the triangle. 


SOLUTION. 
The given triangle is the triangle of reference. The circumconic 


+ vya + wai =o 
will be an equilateral hyperbola if 
ucos A +vcos B+ wcos C=0; 


v7 | - uy u,3 Yva=@ . b . & 


are the equations to the centre; and the resultant of this system is 
+ bya + ca3 = ad cos A + cos B + c7* cos C, 


which represents the nine-points circle. In determinant notation this resultant 


can be written 
O a“ 


| O yr 
icos C cos B cos A oO} 

a form which is probably new. [R. H. Graves. 
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P and Q are middle points of opposite edges of a tetraedron. A plane 
through PQ intersects two other opposite edges in / and .V. Show that JZV 
is bisected by PQ. [Asaph Hall. ] 


SOLUTION I. 


Draw a plane (A) through an edge of the tetraedron and the opposite middle 
point. This plane divides the solid into equal parts, since the parts stand on the 
same base and have equal altitudes. Draw another plane (#4) through the points 
PMQN. This plane cuts off from the equal parts two small tetraedrons which 
are also equal. To show this draw parallel planes through opposite edges 
of the given tetraedron, which will form a circumscribed parallelopipedon, the 
opposite edges of the tetraedron being the cross diagonals of the faces of the 
parallelopipedon. From this construction it follows that the points J7 and V 
are equally distant from the plane (4). Hence the small tetraedrons are equal, 
since they stand on equal bases and have equal altitudes. From this equality it 
follows that the line JZ is bisected by PQ, round which the plane (/) revolves. 

Hall.] 
SOLUTION II. 


Let €,, €, €,, and e, be the vertices of the tetraedron; whence 
P= + &), 
Q= + &), 
M = fie, + (1 — 
and N= ae, + (1 — a)e,. 
M, N, P, Q are coplanar; hence 
(€, + €) + 4m + €,) + n[ae, + (1 — a) = fe, + (1 — &. 
Reducing, we have 
(}/ + na) e, + (4/— + (4m — 1 4+ =[n(a— 1) — e,, 


a numerical relation between four points not coplanar; hence the coefficients of 
the points are each equal to zero; whence a = 3, and 


(M+ N) =$(ae, + +4 [0 —A 
lies on the line PQ. [/. N. James.] 
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Prof. Graves finds the result anticipated in Byerly’s Integral Calculus: 
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Ir a? = & — bc and P = 2 — ca, then 2 = a* — ad. [Frank Morley.] 


230 


SOLVE the equations 
— z)=a, 


=e. [Frank Morley.] 
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Ir we take products of x consecutive terms of the arithmetical series a, 
a — d, etc., commencing for the first product with the first term, for the second 
product with the second term, and so on; and then multiply these products by 
the coefficients in the expansion of (1 — x)", the aggregate will be ~! d", which 
is independent of the first term. [W. W. Johnson.) 
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IF a square be inscribed in the face of a cube, the plane determined by one 
side and the corner of the opposite face corresponding to the adjacent corner of 
the same face touches the inscribed sphere. [Z7. M. Blakslee.| 
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Ir squares be inscribed in the faces of a cube, their sides and the corre- 
sponding vertices of the opposite faces determine 24 planes tangent to the 
inscribed sphere. M. Blakslee.| 


234 
In a regular heptagon ABCDEFG, show that 


I 


AC* AD [Frank Morley.] 
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Finp the areas of the greatest and of the least rhombus inscribed in an 
ellipse. [R. H. Graves.] 
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